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Abstract 



^ ■ We take a critical view at the basic definition of extended single particle states in a non- 

00 . translationally invariant system. For this, we present the case of a hierarchical lattice and 

incorporate long range interactions that are also distributed in a hierarchical fashion. We show 

(— I ' that it is possible to explicitly construct eigenstates with constant amplitudes (normalized to 

^ ■ unity) at every lattice point for special values of the electron-energy. However, the end-to-end 

transmission, corresponding to the above energy of the electron in such a hierarchical system 

depends strongly on a special correlation between the numerical values of the parameters of 

the Hamiltonian. Keeping the energy and the distribution of the amplitudes invariant, one can 

transform the lattice from conducting to insulating simply by tuning the numerical values of the 

long range interaction. The values of these interactions themselves display a fractal character. 
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1 Introduction 

The problem of localization of single particle states that was initially raised and solved by 
Anderson and others [1] still remains very much alive [21 [3] and, have given birth to an enormous 
and highly rich literature [1]-[I3]. The fundamental result is that, in one and two dimensions all 
the single particle eigenstates are exponentially localized in presence of uncorrelated disorder [1] , 
while in three dimensions one witnesses the possibility of a metal-insulator transition. 

Interesting twist in the concept of electron localization in low dimensional systems came 
up with the work of Dunlap et al lU [8] who discussed that a short range positional correlation 
between the constituents of a one dimensional disordered chain could lead to resonant (extended) 
eigenstates with high transmittivity. The basic cause was traced back to a local resonance in 
a cluster of impurities, and the case was referred to as the random dimer model (RDM). The 
idea was tested in low dimensions with various form of disorder by many others [9]-|ll]. Later, 
it was argued, based on numerical diagonalization of the Hamiltonian that, suitable long range 
positional correlations between the potentials assigned to the atomic sites in a linear chain of 
atoms one could initiate a m,etal-insulator transition even in one dimension |121ll3j. 

The idea of generating extended eigenstates in low dimensional systems without any trans- 
lational order (disordered systems provide one such class) was further extended to the infinite 
quasi-periodic chains, where one could unravel an infinite variety of these unscattered states 
owing to the self similarity of the lattices |14t I15| . Similar studies were carried out on fractal 
networks |16j-|20j where, even in the absence of any local resonating clusters such as in the 
case of the RDM [6] or in the case of a quasi-periodic chain [T3] one finds an infinite number 
of extended single particle states. In fact, Schwalm and Moritz [20] have argued, based on an 
extensive numerical support that, even a continuum of extended eigenstates might exist in a 
class of fractal lattices, a possibility that was also pointed out earlier somewhere else |19] . 

In all these exciting conceptual developments, a critical issue is practically overlooked. This 
is the question of categorizing an extended (resonant) state. More explicitly, let us raise the 
question, "when do we call an eigenstate extended ?" . An obvious answer is obtained by looking 
at the amplitude of the wave function at a given energy. A non trivial distribution of the 
amplitude throughout an infinite lattice very legitimately points to an extendedness of the 



eigenfunction. A second way to characterize an extended wave function is through a calculation 
of the electronic transmission across the lattice at the concerned energy. Usually, an extended 
state should give rise to a ballistic transmission. This second criterion is compatible to the first 
one in all systems where translational symmetry prevails. However, its not apparent whether 
these two demands always go hand in hand in systems where extended states exist even in the 
absence of any translational order (such as the fractals for example), and a serious introspection 
of the issue is in order. 

In the present communication, we specifically address this problem. We provide explicit 
example of a hierarchical lattice with long range interactions where, eigenstates with identical 
(and non-zero) amplitudes can be constructed at all the lattice points for a special energy 
of the electron. This construction demands a well defined correlation between the numerical 
values of a subset of the Hamiltonian parameters. The energy of the electron can be chosen 
independent of all or, a sub-set of the hierarchical long range interactions. It is then shown that 
the corner-to- corner propagation of an electron depends crucially on the strengths of the long 
range interactions rendering, in some cases, the lattice completely transparent to an incoming 
electron. In other situations, with a different choice of the hierarchical parameter, a topologically 
identical lattice with the same constant distribution of amplitudes of the eigenfunction, becomes 
completely opaque to an electron with the same energy as in the first case. This observation, to 
our mind, introduces in a possible conceptual confiict between the extendedness of an eigenstate 
and a ballistic transmission, particularly in such hierarchical systems. 

We discuss two kinds of hierarchical interaction in a fractal network designed in the line of the 
well known Berker lattice [2T] . Hierarchical structures using superconducting wire networks have 
already been fabricated and studied experimentally [22]. With the present day nano-technology, 
tailor-made geometries with quantum dots are also possible to fabricate using scanning tunnel 
microscope(STM) as tweezers. Our proposed structures are therefore not far from reality. Also, 
by controlling the proximity of the dots one can induce tunnel hopping almost at will. 

We work within a tight binding approach. A real space renormalization group (RSRG) deci- 
mation scheme [53] is employed to calculate the corner-to-corner electronic transmission in finite 
but arbitrarily large lattices, and to analyze the character of the electron states. Incidentally, 



a similar analysis has been made recently [23] on a Koch fractal with hierarchical interactions, 
and similar questions have been raised. The present lattice is a non-trivial generalization of the 
Koch fractal [25]. and exhibits a richer spectrum of results, yet maintaining the basic issue, viz., 
the nature of extended eigenstates in such fractal lattices. 

In what follows, we describe the results. In section 2, we describe the model and the essential 
method in resolving the problem. Numerical results and related discussion are presented in 
section 3, and we draw our conclusions in section 4. 

2 The model and the method 

2.1 The Hamiltonian and the decimation scheme 

We begin by referring to Fig. [1] and Fig. [2j We have designed a hierarchical lattice with two 
different configurations, in one of which the long range interaction is along the axial direc- 
tion (Fig. [1]) and in the other one, it is along the transverse direction (Fig. [2]). We adopt a 
tight binding formalism, and incorporate only the nearest neighbor hopping. The tight binding 
Hamiltonian for non- interacting electrons in the Wannier basis is given by, 

H = Y^e,\i){i\+Y,[tij\i){j\ +tji\j){i\ ] (1) 

where, e^ is the on-site potential of an electron on the i-th. atomic site and tij = tji is the nearest 
neighbor hopping integral between the i-th and j-th sites. For the nearest neighboring sites 
of the lattice we assume tij = t, while the long range hopping integrals are tij = r, and are 
assumed to follow a particular rule: t(?i) = A"t where, n represents the level of hierarchy, and 
A is the hierarchy parameter. For such a hierarchical lattice of finite but arbitrarily large size, 
the on-site potential e^ will be assigned values e^ for the two extreme sites and e5(„) , ec{n) ^tc. 
for the bulk sites depending on their positions on the lattice as explained in Fig. [1] and Fig. [2l 
Exploiting to the self similarity of the lattice, we can apply the real space renormalization 
group (RSRG) decimation scheme [23] to decimate an appropriate subset of sites. In decimating 
those subset of sites, we have used the standard difference equation, which is an equivalent 



description of the Schrodinger equation, viz., 

j 

Here, ^/Ji is the amphtude of the wave function at the i-th atomic site and the index j represents 
the nearest neighbors of i. The recursion relations for the site energies and hopping integrals 
for the two cases are given by: 

I. The axial case 
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where, S = E — ec(i) ~ ''"(1) ^^i^l fi = 5[E — £^(1)] ~ 2t^ (6) 

We present our numerical results in the next section using the above set of recursion relations. 



2.2 The transmission coefficient 

To get the end-to-end transmission coefficient of an i-th generation fractal, we clamp the system 
between two semi-infinite ordered leads. The leads, in the tight binding model, are described by 
a constant on-site potential eg and a nearest neighbor hopping integral to- We then renormalize 
the lattice i times to reduce it into an effective dimer consisting of the 'renormalized' j4-atoms, 
each having an on-site potential equal to e)^ and the effective A-A hopping integral t^^'. The 
transmission coefficient is then obtained by the well known formula [26], 

^ 4sin^A;a 

T = (7) 

[{Pu - P21) + {Pn - P22) cos kaf + [(Pn + P22) sin kaf 

where, 
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P21 = —P12 

P22 = -j^ (8) 

Here, cos ka = {E — eo)/2to 

a is the lattice spacing in the leads, and throughout the calculation we shall set eg = 0, and 
to = l. 



3 Numerical results and discussion 

3.1 The eigenvalue spectrum 

Before discussing the precise point of interest, we prefer to have a look at the eigenvalue spectra 
of the proposed hierarchical structures. In particular, we examine the formation of the bands 
and the gaps as a function of the hierarchy parameter A both for the axial and the transverse 
cases. To this end, we use a well known trick used frequently in the case of a quasi-periodic 
lattice [23 , and later used for hierarchical structures as well \2A\ [2H] . We sequentially generate 
periodic approximants of the original hierarchical structure, calculate the trace of the transfer 
matrix corresponding to a 'unit cell' and extract those values of energy for which the magnitude 
of the trace of the transfer matrix remains bounded by 2 [27]. The results are shown for the 
axial and the transverse cases in Figs. [3ja) and (b) respectively. The common feature in both 
the figures is the presence of multiple bands and gaps. A variation of the hierarchy parameter A 
leads to band overlapping. Specifically speaking, in the axial case, the density of allowed energy 
values is larger in the range — 0.2t < A < 0.2t. Band crossings maximize in this area. An 
increase in the numerical value of A leads to a thinning of the spectrum. Influence of A on the 
spectrum of the transverse model is also similar. 

3.2 The unusual states 

We now examine some special situation which is the focus of this article. 

(a) The axial model 

Let us begin with the axial model. It is easy to verify that, if we choose the energy of the 
electron E = eA + t, then one can construct, by hand, a wave function with amplitude equal to 
unity at every lattice point, provided one also fixes 

eB(n) +T{n) = EA- 2t 

^C{n) = ^A-t for all n > 1 (9) 



Two points are worth noting. First, the distribution of amphtude thus constructed is inde- 
pendent of the individual values of e^fn) ^ind T{n), and only requires the special correlation 
in their numerical values as suggested in Eq. ^. Thus, in principle, eB{n) can be chosen out 
of any sequence of uncorrelated random numbers. The values of the hopping integrals T{n) 
only then need to be selected accordingly. Thus the constructed amplitude distribution remains 
valid even for a random choice of a subset of the on-site potentials eBtn)^ ^iid presents a new 
kind of extended wave function. However, this construction can not automatically be related to 
the good transmission property of any large but finite lattice. This is easily understood when 
one appreciates that the choice of the energy didn't depend on the hierarchically distributed 
hoppings viz., T{n), but the end-to- end transmission is bound to be sensitive to the individual 
values of r(?i). 

This is illustrated in Fig. Hl^a), where we have worked out the transmission coefficient across 
a 7th generation hierarchical network and for T{n) = \T[n — 1), with t(1) = \t. Energy is set 
aX, E = eA + t. The corresponding eigenstate has equal amplitude (normalized to unity) at all 
lattice points. The transmission spectrum shows that the 7th generation network is transparent 
to an incoming electron with the above energy only when A assumes a specific set of values. The 
transmission in this range of A is completely ballistic for certain values. We may assign the wave 
function in this case the status of an extended state. Interestingly, the same distribution of the 
amplitudes of the wave function, at the same energy E = eA + t makes the lattice completely 
opaque to the incoming electron for other ranges of A. The energy definitely corresponds to an 
eigenstate of the system, as has been verified by calculating the local density of states at the 
sites of an infinite hierarchical lattice for the axial case. Also, the very fact that, one is able to 
construct such a state on a lattice, no matter how large it is, automatically confirms that it is 
an eigenstate. This observation leads to the question of a proper identification of an extended 
state in a non-translationally invariant system. 

The procedure can be implemented on say, a one step renormalized lattice. That is, we can 
extract an energy eigenvalue by solving the equation E = e^ + f'. The energy turns out to be 
a function of ea-, t, ec(i); ^-B(i) ^^"^ ''"(1)) but remains independent of es^n)-, ^C{n) arid r(n) for 



n > 2. Once again we are to choose, 

eB(n) + T{n) = e^ - 2t' 

e^(„) = e'^-t' for all n > 2 (10) 

The entire scheme works as before. In Fig. [5)^ a) we illustrate the transmission coefficient for 
E = 0. This energy is extracted by solving the equation E = e'^ + t' , where we have chosen 
^A = ^B{i) = ^C(i) = Oi ^iid t = r(l) = 1. A fine scan of a selected portion of Fig. [5]J^a) is 
presented in Fig. [5]^b) to show the self similar distribution of the values of A. The process can, 
in principle, be continued and one can extract energy eigenvalues for such unusual eigenstates 
by solving the equation E = e^ + t^ ' at any ^-th stage of renormalization. Ideally, we thus 
have an infinite number of such eigenstates. 

(b) The transverse model 

For the transverse model, we can use an identical string of arguments to construct an eigenstate 
with an amplitude equal to unity at every lattice point. For this construction we fix the energy 
E = eA + t, and demand, 

ec{n) + t{''T') = ^A — t for all n > 1 (11) 

This choice of parameters leads to a completely new scenario in comparison to the axial case. 
Here, for values of the hierarchy parameter A > 1, it is found that at any i-th stage of renor- 
malization the on-site potential at the C-sites grow following the rule: 



where, c^'(^) is a constant, function of A, and the on-site potential at the -B-sites at any 
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th stage of renormalization eJ, ^ remains constant. Thus an incoming electron with energy 



E = €A + i will face effectively higher and higher potential barriers, offered by the 'C sites while 
travelling through the lattice. This will lead to a gradual decay of the end-to-end transmission 
as the system grows in size. In the thermodynamic limit the hierarchical lattice will remain 
non-conducting. For A < 1, no regular pattern in the on-site potentials is observed. However, 
the hopping integral t always decays to zero for E = eA + t. 

We have also examined the case, where the energy is extracted from a one step renormalized 
lattice by solving the equation E = e'^ + t'. For this we are free to choose ca, e_B(i)i fc(i)) 
t and r(l) freely, that is, in an uncorrelated fashion. The correlation now sets in from the 
hierarchy level n = 2 onwards. For example, we have examined the special situation when 
eA = eB(i) = ec(i) = 0, and t = r(l) = 1. The roots are, E = -1.90321, 0.193937, and 2.70928. 
The transmission in each case drops fast as the lattice grows in size. The amplitude-distribution 
on a one step renormalized lattice is such that ipi = 1 on every vertex of the renormalized lattice. 
Therefore, the states are still strictly localized. Thus we can say that, a hierarchically distributed 
long range hopping in the transverse direction does not allow the lattice to be conducting. 

Thus, once we appreciate that a long range tunnel hopping can be associated with the proxim- 
ity of the atoms in the structure, we see that a proximity along the principal axis of the fractal 
allows both for conduction and localization, whereas, a proximity in the transverse direction 
makes the lattice non-conducting in general. 

Before we end, it should be emphasized that, all these discussions are made with reference 
to special values of the energy E for which one can construct a unique spatially extended 
distribution of the eigenfunctions. However, hierarchical lattices, as already discussed in the 
introduction, may possess both localized and extended eigenstates coexisting in the spectrum. 
For example, we have worked out a special situation in the transverse case, where a different set 
of values for the Hamiltonian parameters and the electron-energy may lead to a one cycle fixed 
point of the parameter space. The eigenstate in this particular case is definitely extended and 
the transport is high. However, we refrain from entering into this aspect to save space. 

10 



4 Conclusions 

In conclusion, we have undertaken a detailed study of the electronic states and transport across 
a hierarchical lattice corresponding to a special set of energy eigenvalues. The central issue has 
been to examine if an unambiguous definition of the extendedness of an eigenstate in a lattice 
that lacks translational symmetry is obtainable. We come to the conclusion that, a state that 
looks extended by construction is not necessarily conducting, and the mobility of the state is 
strongly sensitive to a correlated choice of a subset of the system parameters. Even within 
the same basic lattice topology, a long range correlated tunnel hopping along the principal axis 
is found to lead to both extended and localized states depending on the value of the hierarchy 
parameter. The hierarchically long range hopping in the transverse direction, on the other hand, 
makes the eigenfunction localized in the lattice. 
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Figure Captions: 

Figure 1: (Color online), (a) Schematic diagram of the 2nd generation of a hierarchical lattice 
with long range interaction along the 'axial' direction. 'A' denotes the extreme sites, while 
B{n) and C{n) represent the bulk sites depending on their positions in the lattice, n being the 
hierarchy index, 't' represents the nearest neighbor hopping integral and T{n) represents the 
value of the long range hopping integral in the n-th level of hierarchy, (b) The renormalized 
version of (a). 

Figure 2: (Color online), (a) Schematic diagram of the 2nd generation of a hierarchical lattice 
with long range interaction along the 'transverse ' direction. 'A ' denotes the extreme sites, while 
B{n) and C{n) represent the bulk sites depending on their positions in the lattice, n being the 
hierarchy index, 't' represents the nearest neighbor hopping integral and T{n) represents the 
value of the long range hopping integral in the n-th level of hierarchy, (b) The renormalized 
version of (a). 

Figure 3: (Color online). Energy eigenvalue spectrum of a hierarchical fractal lattice as a 
function of the hierarchy parameter A, obtained from the trace of the transfer matrix for the 
7th generation fractal, taken as the 'unit cell'. We have set ea = ^B(n) = ^C(n) = with n = 7 
and t = 1. (a) The axial case and (b) The transverse case. 

Figure 4: (Color online), (a) Transmission coefficient across a 7th generation fractal network 
(the axial case) for £■ = e^ + 1 with eB(n) + ''"("-) = ^A — 2t, and ec(n) = ^A — t. (b) Fine scan of 
a selected part of (a) to reveal the self-similar distribution of A. We have set ca = 0, and t = 1. 

Figure 5: (Color online), (a) Transmission coefficient across a 7th generation fractal network 
(the axial case) for E = e'j^ + t' and (b) a fine scan of a selected part of (a) to reveal the 
self-similar distribution of A. We have set e^ = ^B{i) = ^C(i) = 0) ^'^'^ * = I'iX) — 1- 
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Figures: 
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Figure 1: (a) Schematic diagram of tlie 2nd generation of a hierarcliical lattice with long 
range interaction along the 'axial' direction. 'A' denotes the extreme sites, while B[n) 
and C{n) represent the bulk sites depending on their positions in the lattice, n being the 
hierarchy index, 'i ' represents the nearest neighbor hopping integral and T{n) represents the 
value of the long range hopping integral in the n-th level of hierarchy, (b) The renormalized 
version of (a). 
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Figure 2: (a) Schematic diagram of the 2nd generation of a hierarchical lattice with long 
range interaction along the 'transverse ' direction. 'A 'denotes the extreme sites, while B{n) 
and C{n) represent the bulk sites depending on their positions in the lattice, n being the 
hierarchy index, 't ' represents the nearest neighbor hopping integral and r (ra) represents the 
value of the long range hopping integral in the n-th level of hierarchy, (b) The renormalized 
version of (a). 
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Figure 3: Energy eigenvalue spectrum of a iiierarchical fractal lattice as a function of the 
hierarchy parameter A, obtained from the trace of the transfer matrix for the 7th generation 
fractal, taken as the 'unit cell'. We have set e^ = (^(^n) = ^C{n) = with n = 7 and t = 1. 
(a) The axial case and (b) The transverse case. 
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Figure 4: (a) Transmission coefficient across a 7th generation fractal network (tlie axial 
case) for E = eA + t with eB{n) + ''"('T') = ca — 2t, and ec(n.) = eA — *• (b) Fine scan of a 
selected part of (a) to reveal the self-similar distribution of A. We have set e^i = 0, and 
t = 1. 
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Figure 5: (a) Transmission coefficient across a 7th generation fractal network (the axial 
case) for E = e'^ + t' and (b) a fine scan of a selected part of (a) to reveal the self-similar 
distribution of A. We have set e^ = £5(1) = ^C{i) = Oj ^^^ * = ''"(1) = 1- 
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